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Abstract 


As an important duality result in linear optimization, the Goldman—Tucker 
theorem establishes strict complementarity between a pair of primal and 
dual linear programs. Our study extends this result into the framework of 
linear fractional optimization. Associated with a linear fractional program, 
a dual program can be defined as the dual of the equivalent linear program 
obtained from applying the Charnes—Cooper transformation to the given 
program. Based on this definition, we propose new criteria for primal and 
dual optimality by showing that the primal and dual optimal sets can be 
equivalently modeled as the optimal sets of a pair of primal and dual lin- 
ear programs. Then, we define the concept of strict complementarity and 
establish the existence of at least one, called strict complementary, pair 
of primal and dual optimal solutions such that in every pair of comple- 
mentary variables, exactly one variable is positive and the other is zero. 
We geometrically interpret the strict complementarity in terms of the rel- 
ative interiors of two sets that represent the primal and dual optimal sets 
in higher dimensions. Finally, using this interpretation, we develop two 
approaches for finding a strict complementary solution in linear fractional 
optimization. We illustrate our results with two numerical examples. 
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1 Introduction 


A mathematical optimization problem is specified as a (primal) linear frac- 
tional program (LFP) when a linear fractional function (i.e., ratio of two 
affine functions) is optimized subject to a set of linear constraints on the given 
variables.!_ The linear fractional optimization frequently appears in a wide 
variety of real-world applications, including information theory, numerical 
analysis, game theory, cutting stock problems, shipping schedules, macroeco- 
nomic planning model, and so on. More details can be found in [1, 11, 24, 25] 
and references therein. It is also applied in the measurement of efficiency by 
data envelopment analysis; see, for example, [10, 15, 22, 26, 27, 28, 29, 30, 31] 
among others. Therefore, considerable research interest has been devoted to 
this branch of optimization. 

The literature on the duality of linear fractional optimization associate 
various duals to the primal LFP. Chadha [7] suggested a dual in the form of a 
linear program (LP) and proved some duality statements directly. However, 
the constant scalars in the numerator and denominator of the primal objective 
function are assumed to be absent in their work. Chadha and Chadha [8] 
extended Chadha’s results to the general case, where the constant scalars 
are taken into account. An interesting note regarding their impressive work 
is that their results can be deduced in an alternative way from the duality 
of linear optimization. In fact, an indirect approach for constructing a dual 
program is to transform the primal LFP into an equivalent problem that its 
dual can be constructed in the classical way; see [25]. By using this approach, 
it can be verified (as in Section 2) that the dual program proposed in [8] is 
nothing else than the dual of the equivalent LP resulting from applying the 
well-known transformation of Charnes and Cooper [9] to the primal LFP. 

Though demonstrating the common complementary slackness condition 
between the primal LFP and its dual, Chadha and Chadha [8] did not inves- 
tigate the strict complementarity between them. Furthermore, to the best 
of our knowledge, no other research exists on such investigation. Motivated 
by these, we extend an important duality result proved by Goldman and 
Tucker [12] from linear optimization to linear fractional optimization. The 
so-called Goldman—Tucker theorem establishes the strict complementarity 
between a pair of primal and dual LPs. It states that at least one, so-called 
strict complementary, pair of primal and dual optimal solutions exists such 
that the sum of each pair of complementary variables is positive. That is, in 
every pair of complementary variables, exactly one variable is positive and 
the other is zero; see, for example, [20] for more details on the theory and 
applications of strict complementarity in linear optimization. 


1 Tf the given objective function is optimized with no restrictions on the values of its 
variables, then the optimization problem is called unconstrained. Useful information on 
approaches developed for solving unconstrained optimization problems can be found in 
[3, 16, 21], among others. 
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As a complementary to the work of Chadha and Chadha [8], this paper 
shows that the primal and dual optimal sets can be equivalently modeled as 
the optimal sets of a pair of primal and dual LPs. Using this fact, we propose 
new criteria for primal and dual optimality in terms of the belongingness of 
these LPs’ objective vectors to the binding polyhedral cones at primal and 
dual feasible solutions. Then we define the strict complementary slackness 
condition for an LFP and demonstrate the existence of a strict complemen- 
tary solution. We also show that any strict complementary solution induces 
unique optimal partitions for the sets of indices of nonnegative decision vari- 
ables. 

To deal with the problem of finding a strict complementary solution, we 
equivalently represent the primal and dual optimal sets by two nonnegative 
polyhedral sets in higher dimensions, which are described only by equality 
defining constraints.? Then we geometrically interpret the strict complemen- 
tarity by proving that any pair of relative interior points of these polyhedral 
sets is a strict complementary solution, and vice versa. Based on this inter- 
pretation, we turn the problem under consideration to the equivalent problem 
of identifying a maximal element of a nonnegative polyhedral set. Exploiting 
the recent work of Mehdiloozad et al. [20], who have addressed the latter 
problem, we develop two linear optimization approaches for finding a strict 
complementary solution. 

The remainder of this paper is organized as follows. Section 2 provides 
the necessary background needed for the rest of the paper. Section 3 pro- 
poses new criteria for primal and dual optimality and illustrates them with 
a numerical example. Section 4 establishes the strict complementarity for 
LFPs. Section 5 proposes an LP for finding a maximal element of a non- 
negative polyhedral set and, thereby, develops two approaches for finding a 
strict complementary solution. Section 6 illustrates these approaches by a 
numerical example. Section 7 contains concluding remarks and suggestions 
for future research. Appendix A provides the GAMS (General Algebraic 
Modeling System) code of our proposed approaches. 


2 Background 


2.1 Notation 


Let R¢ denote the d-dimensional Euclidean space, and let Re denote its non- 
negative orthant. We denote sets by uppercase calligraphic letters, vectors by 
boldface lowercase letters, and matrices by boldface uppercase letters. We 


2 A polyhedral set is said to be nonnegative if it is a subset of the nonnegative orthant 
of Euclidean space. By the “defining constraints” of such a polyhedral set, we refer to 
the constraints imposed other than the nonnegativity conditions. 
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denote the cardinality of a set S by Card(S). By convention, all vectors 
are column vectors. The superscript T denotes the transpose of a vector or 
matrix. 


Vectors 0 and 1 are vectors all components of that are equal to 0 and 1, 
respectively. The dimensions of these vectors are clear from the context in 
which they are used. For simplicity, the notation (a;b) € R¢+’ is used to 
show the column vector obtained by adding vector b € R“ below the vector 
a € R¢. For vectors a,b € R%, the inequality a > b (resp., a > b) means 
that a; > b; (resp., a; > b;) for alli =1,...,d. 

Matrix O is the matrix all components of that are equal to 0, and matrix 
I is the identity matrix. The dimensions of these matrices are clear from the 
context in which they are used. We denote the ith (i = 1,...,d) row and the 
jth (j = 1,...,d’) column of a dx d’ matrix A by a’ and aj, respectively. In 
particular, we use the notation e; to denote the jth column of the identity 
matrix of size d x d, that is, e; = (Use, 1 een: € R¢ for j = 1,...,d. 

jth 

Recall from [23] that the relative interior of a subset Y of R?%, denoted by 
ri(%), is defined as the interior we get when 4 is regarded as a subset of its 
affine hull, denoted by aff (4). Formally, 


ri(¥) = {x° € ¥: Nz (x°) Naff (¥) C ¥ for some € > 0}, 


where NV, (x°) = {x € R#: ||x — x°|| < e}. 

Recall also from [20] that any convex (and, in particular, polyhedral) sub- 
set of R4 is called a nonnegative convex (polyhedral) set. Additionally, any 
element of a nonnegative convex set is said to be maximal, if the number of 
its positive components is maximum. We denote the support of a nonnega- 
tive vector a € R4 by supp (a), that is, supp (a) = {i € {1,...,d}: a; > 0}. 
We also denote by me(4) the set of all maximal elements of 1, that is, 


me (4’) = argmax Card (supp (x)). 
xEer 


2.2 Linear fractional program 


A function of variables is said to be linear fractional if both its numerator 
and denominator are affine functions of the given variables. A mathematical 
optimization problem that optimizes a linear fractional objective function 
subject to a set of linear constraints is called as a linear fractional program 
(LFP). Formally, the general form? of the primal LFP is defined as 


3 The standard form of the primal LFP results from (1) by replacing the inequality sign 
“<” in (1b) by the equality sign; see [1, Section 1.3]. 
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c x+a 


max f (x)= ax4B (la) 
subject to 

Ax <b, (1b) 
x>0, (1c) 


where x € R” is the vector of decision variables, c € R” and d € R” are, 
respectively, numerator and denominator vectors of objective function, a € R 
and 6 € R are objective scalars, A is an mx 7 constraint matrix, and b € R™ 
is the right-hand side vector. 


Let S denote the feasible set of program (1), which is clearly a polyhedral 
set in R”. To ensure that the function f is well-defined on S, it is assumed 
that its denominator maintains a constant sign on S. Without loss of gener- 
ality, we assume that d'x+ 6 > 0 for all x € S. Then the objective function 
f is both quasi-convex and quasi-concave over S and, therefore, every local 
maximum is a global maximum (see, e.g., [3]). To guarantee the occurrence 
of finite optimality for program (1), we also assume that S is regular (i.e., 
nonempty and bounded). 


An effective approach for solving program (1) is to transform it into an 
equivalent LP by the well-known Charnes—Cooper transformation [9]. In fact, 
if we define t = ax and x = tx, then multiplying both sides of (1b) by t 
converts program (1) to the following LP: 


max c'X+at (2a) 
subject to 

Ax — bt <0, (2b) 
d'x+6t=1, (2c) 
x>0,t>0. (2d) 


Let S be the feasible set of program (2). Because t > 0 for all (x;t) € S,* 
the following implication between the feasible solutions of programs (1) and 
(2) is established: 


= 1 
(xthEeS = FX ES,t> 0. (3) 


Especially, if (x*;¢*) is an optimal solution to program (2), then ax* is an 
optimal solution to program (1) (see, e.g., [1, p. 57]). 


4 Indeed, if t = 0 for some (x;t) € S, then it follows from (2b)—(2c) that AX < 0 and 
x #0. This means that the vector x is a recession direction of the feasible set S, thereby 
contradicting the regularity assumption. Therefore, t > 0 for all (x;t) € S. 
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2.3 Dual of linear fractional program 
To state the dual to LFP (1), let the vector y € R™ be dual to (2b) and the 


scalar z be dual to (2c). Then, by the duality of linear optimization, the dual 
to LP (2) is stated as follows:° 


min g(y;z)=2 (da) 

subject to 

Aly+dz>c, (4b) 
—b'y + 6z=a, (4c) 

y > 0, z sign free. (4d) 


Observe that program (4) is nothing else than the LP introduced in [8] 
as the dual of program (1). We denote by D the feasible set of this program. 

Throughout this paper, LP (4) is defined to be the dual of LFP (1). 
The next three theorems demonstrate the duality relationships between pro- 
grams (1) and (4). 


Theorem 1 (Weak duality). [8] For any x € S and any (y;z) € D, we have 
f(x) Sg (y;2). 


Theorem 2 (Optimality criterion). [8] If the feasible solutions x € S and 
(y;z) € D satisfy f(x) = g(y;z), then they are optimal solutions to pro- 
grams (1) and (4), respectively. 


Theorem 3 (Strong duality). [8] If x* is an optimal solution to program (1), 
then there exists some optimal solution (y*;z*) to program (4) such that 
f(x") =g(y"32"). 


The following result gives a necessary and sufficient optimality condition, 
called complementary slackness condition (CSC), in terms of the complemen- 
tarity of the primal and dual feasible solutions. 


Theorem 4. [8] Feasible solutions x* € S and (y*;z*) € D are optimal if 
and only if they fulfill the following conditions: 


vil x* = u*ly* = 0, (5) 


where u* = b — Ax* and v* = Aly* + dz* —c. 


From Theorem 4, the pairs (23,05), j = 1,...,n, and (ui, yi), i= 
1,...,m, are called complementary variables. 


5 Note that the inequality constraint —b' y +z > a has been replaced with its equality 
form in program (4), because the optimal value of its corresponding dual variable in 
program (2), t, is always positive. 
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3 Our criteria for primal and dual optimality 


In this section, we propose new criteria for the optimality of LFP (1) and its 
dual (4), and present their geometrical interpretations. 

Denote by f* the optimal objective value of program (1), and let d* = f*d 
and §8* = f*@. Then the set of all optimal solutions of program (1) can be 
defined by conditions (1b)—(1c) and the additional equality requiring that the 
objective function of program (1) to be equal to f*. Equivalently, this set is 
stated as follows: 


ae {x € R”: (c—d*)'x=-a@+ £", Ax <b, x> 0}. 
Similarly, an equivalent statement of the optimal set of program (4) is 
{(yi2") eR™?: b'y=—-at+f*, A'y>e-d*,y> o}. 


where z* denotes the optimal objective value of program (4) and is equal to 
f*. Observe that the last components of all optimal solutions of program (4) 
are equal to z*. Therefore, without losing anything, we can remove the last 
dimension of the optimal set of program (4) by projecting it onto the space 
of y-variables. This results the following set: 


yt = {y ER™: by =-a+ 8", ATy >c—a", y 20}, 

which will be loosely referred to as the optimal set of program (4). 
It is clear that the nonempty optimal sets ¥* and Y* are polyhedral 
subsets of R' and R%’, respectively. The next result shows that these sets 


are interestingly the optimal sets of the LP 


max (c—d*)'x 


subject. to (6) 
(1b) — (1c), 

and its dual 
min bly (7a) 
subject to 
Aly = c— d*, (7b) 
y > 0. (7c) 


Theorem 5. Let #5 and F7, be the optimal sets of programs (6) and (7), 
respectively. Then, #5 = V* and Ff = Y*. 


Proof. Let x € X* and y € Y*. Then X and y are, respectively, feasible 
solutions to LPs (6) and (7) such that (c — d*)' x = by. By the optimality 
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criterion theorem of linear optimization, it follows that x € Ff and y € Fp. 
Therefore, 1* C Fs and Y* C FF. 

Conversely, let X € #5 and y € F4. By the strong duality theorem of 
linear optimization, we have (c — d*)' % = b'¥. Because X* # @ and X* 
F5, the weak duality theorem of linear optimization implies that —a + 6* 
b'y. Similarly, it follows from Y* #4 @ and Y* C F% that (ec — dd)" 
—a+ B*. Consequently, we have (c — d*)' x=b!'y¥ =-—a+ 6*. Therefore, 
x € &* and y € Y*, which, respectively, imply Fp C X* and FA C Y*. 


IN IAIN 


Associated with any feasible solution of an LP, the binding cone is defined 
as the convex cone generated by the gradients of all constraints that are 
binding (active) at that solution. Recall from linear optimization that a 
feasible solution to an LP is optimal if and only if its corresponding binding 
cone includes the gradient of the objective function. Based on this fact, we 
apply Theorem 5 to provide necessary and sufficient geometrical conditions 
for feasible solutions of LFP (1) and its dual (4) to be optimal. 

Let x* be a feasible solution to LFP (1), and let G3 denote the union of 
the gradients of all binding constraints at x*, that is, 

Gp = {(a’)" € R”: a’x* = bi} U{-e; €R": r= O}. 

Furthermore, denote by 6% the binding polyhedral cone generated by G5, 
that is, Bj = cone(Gj), where the operator “cone” denotes the conical hull. 
Because the feasible regions of programs (1) and (6) are equal, x* is a feasible 
solution of LP (6). Therefore, the next corollary follows immediately from 
Theorem 5. 


Corollary 1. Let x* ¢ S. Then, x* € A* if and only if ce — d* € By. 


Similarly, let (y*; z*) be a feasible solution to program (4). Additionally, 
assume that Bt, = cone(G7,), where 


Gi = {aj ER”: a} y* =C; — di} U {e; €¢ R™: yF = 0}. 
Then we obtain the following corollary as a consequence of Theorem 5. 
Corollary 2. Let (y*;z*) € D. Then, y* € Y* if and only if b € By. 


We now present a numerical example verifying Corollaries 1 and 2. 


Example 1. Consider the following LFP: 
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Dy 4 
2 frm) =5 


9 
feu) =5 


5 
f@y%2) =| 


6 
fOum)=¢ 


Figure 1: Multiple optimal solutions of program (8) 


621 + 3%2 + 6 


max f (%1,22) = 5a, + 2x0 +5 (8a) 

subject to 

271 + wo <6, (8b) 
—2%7, +42 < 2, (8c) 

21,22 = 0. (8d) 


A graphical approach for finding optimal solution(s) of two-dimensional 
LFPs is to rotate the level-line around its focus point in positive direction 
(i.e., counterclockwise).° Figure 1 illustrates an application of this approach 
to program (8). The feasible region of program (8) in two dimensions 21 
and 22 is the bounded polyhedral set OABC (shaded in gray), and the focus 
point is F = (—1,0). Therefore, the optimal objective value of program (8) 
is f* = 4. Additionally, the set of all optimal solutions to program (8) is 
the segment AB, which is stated below as all convex combinations of the two 
extreme points A and B of the feasible region: 


X* = {x € R?: (ct, 23) = A(0,2) + (1— A) (1,4), A € [0,1]}. 


Table 1 presents the geometrical investigation of the proposed condition 
of primal optimality in Corollary 1 at four extreme points O, A, B, and 
C,, and one nonextreme point P of the feasible region of program (8). As 
expected, the condition holds at the optimal points A, B, and P, but not at 
the nonoptimal points O and C. 

The dual to program (8) is the following LP: 


6 More details on graphical solution of LFPs involving only two variables can be found 
in [1, Chapter 3]. 
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Table 1: The proposed criterion of primal optimality for program (8) 
x f(x) ce-f(x)d Gp Belongingness of c¢ — f (x) d to Bp 
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min z (9a) 
subject to 
2y1 — 2yo2 + 5z > 6, (9b) 
Yt Yo+2z2 3, (9c) 
—6y1 — 2yo + 5z = 6, (9d) 
Y1, 42 > 0, z sign free. (9e) 


By Theorem 3, the optimal objective value of the dual program (9) is 
equal to that of the primal LFP (8), so 2* = q. By Theorem 4, yf = 0 for 
any optimal solution (y{, y3) of (9) because the point A is an optimal solution 
to program (8) for which the inequality constraint (8b) is strict. Additionally, 
both constraints (9b) and (9c) must be binding at optimality because point 
B with both positive components is an optimal solution to program (8). 
Taking these into account, it follows from the constraints of program (9) 
that y3 = 3 for any optimal solution. Therefore, (yi. y3,2*) = (0, 3, 3) is 
the unique optimal solution of LP (9). The stated facts are observable from 
Figure 2, which draws the feasible region of program (9) in three dimensions 
Y1, Y2, and z as the section LM NK of the two-dimensional hyperplane H = 
{(y;z) € R®: — 6y, — 2y2 +5z = 6}. 


Figure 2: Unique optimal solution of program (9) 


Observe that projecting the unique optimal solution of program (4) onto 
the space of y-variables follows that Y* = { (0, ap As an illustration to 
Theorem 5, Figure 3 shows the singleton set Y* to be the optimal extreme 
point of the following LP: 
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min 6y; + 2y2 (10a) 
subject to 
2y1 — yo > > (10b) 
1 
Yit y2= 3° (10c) 
Yi, yo = 0. (10d) 


Figure 3: Representing Y* as the unique optimal solution of program (10) 


Table 2 geometrically investigates the proposed condition of dual opti- 
mality in Corollary 2 at the two extreme points M and N of the feasible 
region of program (9). While the condition is met by the optimal point N, 
it is not true at the nonoptimal point M. This verifies that the projection of 
the optimal point N onto the space of y-variables is in Y*. 


Remark 1. Corollary 1 suggests a geometrical approach for finding optimal 
solution(s) of program (1). It states that the optimality of any feasible point 
in S is equivalent to satisfying the condition given in Corollary 1. Therefore, 
taking into account the fact that the finite optimum must occur at some 
extreme points of S, optimal solution(s) of program (1) can be found by 
examining the proposed condition only at extreme points of S. (A similar 
approach for finding optimal solution(s) of program (4) can be devised based 
on Corollary 2.) It is important to note that this graphical approach does 
not require the enumeration of all extreme points of the feasible region. 
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Table 2: The proposed criterion of primal optimality for program (8) 
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4 Strict complementarity 


4.1 Strict complementary solution 


By Theorem 4, the CSC requires only that the product of each pair of com- 
plementary variables is zero at optimality. Therefore, not only one of the 
complementary variables must be zero at optimality, but also both are al- 
lowed to take simultaneously zero optimal values. It means that the CSC 
does not imply the positivity of pairwise sum of the complementary vari- 
ables. If such positivity holds for a pair of optimal solutions for the primal 
LFP and its dual, then in every pair of complementary variables, exactly 
one variable is positive and the other is zero. Calling this property as strict 
complementarity, we present the following definition. 


Definition 1. Feasible solutions x** € S and (y**; 2**) € D satisfy the strict 
complementary slackness condition (SCSC), if they fulfill the following con- 
ditions in addition to the conditions given in (5): 


xftve>0, ue +y**>0. (11) 


It is clear that feasible solutions to programs (1) and (4) that satisfy the 
SCSC are optimal. We refer to such a pair of solutions as a strict comple- 
mentary solution and denote it by (x**,y**). By the next result, we prove 
the existence of such a strict complementary solution. 


Theorem 6. The following statements are true: 
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(i) Any strict complementary solution to LPs (6) and (7) is a strict com- 
plementary solution to programs (1) and (4). 


(ii) There exists at least one strict complementary solution to LFP (1) and 
its dual (4). 


Proof. Part (i) Let (x**,y**) be a strict complementary solution to LPs (6) 
and (7). Then, by the definition of strict complementary in linear optimiza- 
tion, x** and y** are, respectively, optimal solutions to these LPs, such that 


yisl ys — ay = 0, x*5 +4 ves 0, u*s +p yy > 0, (12) 


where u*® and v** are, respectively, the slack vectors added to the inequality 
constraints in Y* and y*. 

By Theorem 5, x** and (y**; 2**) are optimal solutions to programs (1) 
and (4), respectively. Furthermore, it follows from (12) that these solutions 
meet the SCSC in the sense of Definition 1. Therefore, (x**,y**) is a strict 
complementary solution to programs (1) and (4). 


Part (ti) Because a finite optimum occurs for LPs (6) and (7), the Goldman- 
Tucker theorem implies the existence of a strict complementary solution to 
these LPs, which is a strict complementary solution to programs (1) and (4) 
by part (i) of the theorem. 


4.2 Optimal partitions 


Let (x**, y**) be a strict complementary solution to programs (1) and (4). 
Then, the supports of vectors x** and v** are disjoint and their union is 
equal to the index set {1,...,n}. Similarly, the supports of vectors y** and 
u** form a partition for the index set {1,...,m}. Formally, we can write 


supp (x**) Msupp(v*?) = @, supp (x**) Usupp (v**) = {1,...,n}; 


*S *S *S *S (13) 
supp (u**) M supp (y**) = @, supp (u**) Usupp (y**) = {1,...,m}. 


We call the above partitions as the optimal partitions induced for pro- 
grams (1) and (4). By the next result, we show that these partitions, being 
independent from the given strict complementary solution, are unique. 


Theorem 7. The optimal partitions induced for programs (1) and (4) as in 
(13) are the same across all strict complementary solutions and are, therefore, 
unique. 


Proof. By contradiction, let (x',y') and (x?,y”) be two distinct strict com- 
plementary solutions such that xs > 0 and ue = 0 for some j € {1,...,n}. 
Then, vs = 0 and v3 > 0. Because both optimal sets of programs (1) and (4) 
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are convex, 5(x!,y') + 4(x?,y”) must be an optimal solution such that 


(ut + u’) and $(v" + v’) are, respectively, its primal and dual slack 
vectors. For this solution, we have the contradiction (with the CSC) that 


5 (a + x?) > 0 and 3 (v3 + v2) > 0. Therefore, the optimal partition of 
{1,...,n} is unique. The uniqueness of the optimal partition of {1,...,m} 


follows from a similar argument. 


It is worth noting that the optimal partitions can be useful in situations, 
where knowing the positivity of a variable in some optimal solution of an LFP 
is concerned with. For example, while the slack-based measure (SBM) model 
of Tone [26] is used for the measurement of efficiency in the field of data 
envelopment analysis, the global reference set (peer group) of an inefficient 
decision making unit can be identified by the optimal partition of the index 
set of intensity vector.’ 


4.3 Geometrical interpretation 


We begin this section by recalling the following definition from [4]. 
Definition 2. Let SC R@. A subset St of R%* is a representing set for 
S, if its projection onto the space of x-variables is exactly S, that is, x € S 
if and only if there exists some s € R@ such that (x;s) € S*: 


S= {xeR?: (x;s) € S* for some sR}. 


By adding slack vectors to the inequality constraints of 7* and )*, we 
define the following nonnegative polyhedral sets: 


att = {(x;u) € R°t™: (c—d*)'x =-a+ 6*, Ax+u=b,x>0, u>0o}, 


yt = f(yiv) eR™: b'y =-a+ B, ATty-v=c-d*,y>0,v>0}. 


By Definition 2, V¥** and Y** are polyhedral representing sets for ¥* 
and )Y*, respectively. The next result shows that projecting the relative 
interiors of these representing sets gives the strict complementary solutions 
of programs (1) and (4). 


Theorem 8. It follows that (x**, y**) is a strict complementary solution to 
LFP (1) and its dual (4) if and only if (x**; u**) € ri(4**) and (y**;v**) € 
ri(Y**). 


7 For more details on the concept of global reference set and its identification, the reader 
may refer to [17, 18, 19]. 
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Proof. Let (x**;u**) € ri(¥**) and let (y**;v**) € ri(J**). By [20, The- 
orem 4.1 ], the relative interior of a nonnegative polyhedral set with equal- 
ity defining constraints consists of its maximal elements. It follows that 
(x*8; u*s) € me (A’**) and (y*8;v*8) € me()**). By the Goldman—Tucker 
theorem, (x**, y**) is thus a strict complementary solution to LPs (6) and (7). 
By part (i) of Theorem 6, it follows that (x**, y**) is a strict complementary 
solution to programs (1) and (4). 

Conversely, let (x**, y**) be a strict complementary solution to pro- 
grams (1) and (4). Then, similar to the proof of Theorem 6, it can be proved 
that (x**, y**) is a strict complementary solution to LPs (2) and (4). 


To illustrate the concept of strict complementarity, we return back to 
Example 1. Adding the nonnegative slack variables u,; and uz to (8b) and 
(8c) obtains the following representing set for ¥*: 


Xt = {(x*;u’) ER4: (27,23, up, u2) = (1— A, 4 — 2A, 4,0),  € [0, 1]}. 


Similarly, adding the nonnegative slack variables v, and v2 to (10b) and 
(10c) results the following representing set for Y*: 


Sa {(y:v) ER: (y1,Y2, 01,02) = (0, ;0.0)}. 


Consider the midpoint P = (5,3) of the line segment AB in Figure 1. 


This point is associated with the vector (x? ; uw?) a ((3. 3)" (2, 0)"), which 
is a maximal element and, therefore, a relative interior point of the set **. 
Furthermore, consider the point (0, z) in Figure 3 that is associated with the 
single element (y;v) = ((0. 4)" ;0) of the set Y**. Clearly, (x?,y) is a 
strict complementary solution to programs (8) and (9). 

Note that (x*,y) is a strict complementary solution for all A € (0,1). 
However, this is not true for \ = 0,1. This is because in either of these two 
cases, pairwise sum of the complementary variables is not positive. 


5 Finding a strict complementary solution 


5.1 Finding a maximal element of a nonnegative 
polyhedral set 


Consider the following nonempty polyhedral set in R?: 


P= {xe R*: Px+ Qy+Rz=t, x,y >0, zsign free} , 
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where x € R¢, y € R°, and z € RS are the vectors of variables, P, Q, and R 
are, respectively, matrices of coefficients of orders c x d, cx e, and cx f, and 
t € R° is a constant vector. 


Mehdiloozad et al. [20] developed a general convex optimization program 
for finding a maximal element of a nonnegative convex set. As a consequence 
of their Theorem 3.2, the following result develops an LP for finding a max- 
imal element of P. 


Theorem 9. Let (x!*,x?*, y*,z*,w*) be an optimal solution to the following 
LP: 


max 1!'x! 

subject to 

P (x' +x’) + Qy + Rz=tu, 
1>x'>0, x’,y >0, zsign free, w > 1. 


Then + (x!* + x**) € me(P). 


Proof. By [20, Definition 2.5], the characteristic cone of the nonnegative poly- 
hedral set P is Cp = {ear G) >x € Py, Lay1 > o}. To find a maximal 


element of P, this cone is incorporated into the convex program proposed in 
(20, Theorem 3.2]. This leads to the following LP: 


max 1!x'+w! 

subject to 

P (x! + x’) +Qy+Rz=t(w'+w’), 
1>x'>0,x’,y > 0, z sign free, 1>w'>0, w? >0. 


Because the maximization linear program (16) is feasible and its objective 
function is upper bounded by d+ 1, it has a finite optimal solution, namely, 
a, x?* y* : z*, wh, w*) . 


By the assumption, we have P #4 @. Hence, it follows from [20, Theo- 
rem 3.2] that w'* = 1. It is clear that program (15) is derived from pro- 
gram (16) by replacing w! with its optimal value and using the variable 
substitution w = w! +1. This implies that any optimal solution of pro- 
gram (15) gives an optimal solution to program (16). Namely, if we define 
xl = xl) x?! — x oy’ = y, 2 = 2, wl! = 1, and wu” = w* — 1, then 
(xt, x?’, y’,2’,w'’, w?’) is an optimal solution to program (16). Therefore, 
the statement of the theorem follows from [20, Theorem 3.2 ]. 
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5.2 Our proposed approaches 

Though Theorem 6 demonstrates the existence of a strict complementary 
solution for programs (1) and (4), it does not specify how to identify such a 
solution. To deal with this issue, we develop two approaches in this section 


by applying Theorem 9. The GAMS code of these approaches is provided in 
Appendix A. 


5.2.1 First approach 


From Theorem 8, any pair of relative interior points of Y** and Y** deter- 
mines a strict complementary solution to programs (1) and (4). By Theo- 
rem 9, we develop the following two LPs to find such relative interior points: 


max 1!x! + iu! 


subject to 
(c—d*)" OT] (xi +22) _ (-a+ 6" 17 
| A I u 7 Ww — b WP, ( ) 


x! x? 
> > > ae 
1 = @ = 0, uw ze 0, WP a 1 


max ty spy 


subject to 
b! o! Toy 2 =o * 
K Z| (%, Hee ~ eae wae ") 
1 2 
1> (1) > 0, (2) >0, wp 21. 


Let (x’*,x™*,u’*,u?*,w) and (y'*,y**,y*, y?*, wy) be optimal solu- 
tions to programs (17) and (18), respectively. Then, it follows from Theo- 
rem 9 that 


a= (x!* +x™;ul* +u™) €me(4**), oe) 


(y"sv") = — (y* +y™;vl* +v*) €me("*). (19b) 


* 
Wp 


Therefore, taking into account [20, Theorem 4.1] and Theorem 8, it follows 
from (19) that (x",y") is a strict complementary solution to programs (1) 
and (4). 
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5.2.2 Second approach 


Our first approach of identifying a strict complementary solution requires the 
knowledge of the optimal objective value of program (1). In this section, we 
propose an alternative approach that is exempt from this requirement. We 
exploit the fact that the optimality of feasible solutions to a pair of primal and 
dual LPs follows from the equality of their corresponding objective function 
values. Specifically, we consider the following set: 


W* = {(&y) ER": c K+ at =z, 
Ax < bt, d'x+ 6t=1, 
—b'y+ 6z=a, Aly+dz>c, 
x>0,t>0,y => 0, z sign free}. 


(20) 


By the projection lemma,® it follows that W* is a nonnegative polyhedral 
set in Roe, By adding slack vectors @ and v to the inequality constraints 
of this set, we obtain the following set: 


Were {(%y; tj;v) € R2™+™) : ¢'x + at —z=0, 
Ax+ua—bt=0,d'x+ft=1, 
—b'y+6z=a, A'y+dz-—v=c, 


x,v>0,y,u>0,t> 0, z sign free}. 
(21) 


By Definition 2, W** is a polyhedral representing set for W*. Let 
(X;y;U;v) € W*t. Then (X; y; 0; v) satisfies (21) with some scalars t and z. 
The set W* is defined by conditions (2b)—(2d) and (4b)—(4d) and the addi- 
tional equality requiring that the objective function of LP (2) to be equal to 
the objective function of its dual (4). By the optimality criterion theorem of 
linear optimization, it follows that (x,t) and (y, z) are optimal solutions to 
LPs (2) and (4), respectively. Additionally, u and v are their corresponding 
slack vectors added to inequalities (2b) and (4b), respectively. Consequently, 
we have }(x;a) € X** and (y;v) € Y**. Based on this, the next result 
shows that any maximal element of W*+t determines two relative interior 
points of Y** and Y**, and therefore a strict complementary solution to 
programs (1) and (4). 


Theorem 10. Let (x™°; y™°; w™°; v™°) € me (W**) satisfy (21) with some 
scalars #™° and z™°. Then ( a xe, yo) is a strict complementary solution 
to programs (1) and (4). 


8 The projection lemma states that the projection of a polyhedral set onto the space 
of any subset of its characterizing variables is a polyhedral set; see, for example, [6, 
Corollary 2.4]. 
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Proof. Let (x™°; y™°; w™°; v™°) be a maximal element of W** that satis- 
fies (21) with some scalars t™° and z™°. Then zks (x™°;a@™°) € me (A’**) 
and (y™*; v™°) € me()**). Therefore, by [20, Theorem 4.1], the statement 
follows from Theorem 8. 


Based on Theorem 9, we develop the following LP to find a maximal 
element of W*t: 


max 1s a ea gy ey 


subject to 
c'0 0 0 af ait a —1 0 
x +X 
AI 0 0//( 477) —b 0 0 
d' 0 0 0O a ee ea ed ce ae 
0 0-b' oO} |v, 0%, 0 B a (22) 
A'O 0 -I 0 d c 
ul Ww 
1> vi > 0, xy >0, t>0, z sign free, w > 1. 
vi v? 


Let (aor x a, a*, yl’ 


to program (22), and define 


Ly, v'*,v**,t*, 2*,w*) be an optimal solution 


I 
(Py ay) _— f ( 4 tea 4 you Eile u*; vi* 4 v*) : 
By Theorem 9, we have (x™°;y™°;u'™*;v™°) © me(W*T). If 
ime = +, then it follows by Theorem 10 that (;2x™°,y™) = 


(4 (x +x"), 4 (y* +y*)) is a strict complementary solution to pro- 


grams (1) and (4). 


6 Numerical example 


In this section, we illustrate our proposed approaches of finding strict com- 
plementary solutions with a numerical example, taken from [1, 14]. 


Example 2. Consider the following LFP: 
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%1+2%9+3.543+ w4+1 
221 + 272 + 3.543 + 844 +4 
subject to 


max 


24, + %24+3%34+ 3x4 < 10, (23) 
1, +2%2+ 3+ a4 < 14, 
©1,02,%3,24 = 0. 
The dual of program (23) is the following LP: 
min z 
subject to 
2yi+ yet 2z>1, 
+ Qyo+ 22> 2, 
Y1 Y2 2 (24) 


38yi+ yo+3.52z > 3.5, 
841+ Yyot 3221, 
—10y, — 14y2 + 42 = 1, 
Y1,y2 = 0, z sign free. 


Let us add primal slack variables u; and ug to program (23), and dual 
slack variables v1, v2, v3, and v4 to the inequality constraints of program (24), 
to turn their inequality constraints to equalities. We use our proposed ap- 
proaches to find a strict complementary solution to the above pair of pro- 
grams. Running the modified version of the GAMS code provided in Ap- 
pendix A results that the joint optimal objective value of programs (23) and 
(24) is equal to 0.857. Furthermore, the strict complementary solution ob- 
tained from both of our proposed approaches is 


x; =0, v5 =1.071, 73 =1.2, 23 =0, ui = 0, us = 0; 


vt =1.071 vs =0 vs =0, vx =2.071, ys =0.143 yx = 0.071. 


7 Concluding remarks 


An indirect approach for establishing duality results in linear fractional op- 
timization is based on applying the well-known transformation of Charnes 
and Cooper [9]. This approach converts a primal LFP into an equivalent 
LP and then defines the dual of the obtained LP as the dual of the primal 
LFP. An advantage of using this approach is that it allows exploiting the 
duality results of linear optimization for establishing duality statements in 
linear fractional optimization. 

In this paper, we show that the dual program derived from the above 
approach is the same dual program suggested in [8]. Based on this version of 
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duality, we provide new criteria for primal and dual optimality as our first 
contribution. We equivalently represent the primal and dual optimal sets 
as the optimal sets of a pair of primal and dual LPs. By this representa- 
tion, it follows that a primal (resp., dual) feasible solution is optimal if and 
only if its binding polyhedral cone contains the objective vector of the corre- 
sponding primal (resp., dual) LP. This condition not only is (theoretically) 
necessary and sufficient for the optimality of any general LFP, but also is a 
new geometrical tool for solving two- and three-dimensional LFPs. 

As our second contribution, we introduce the concept of strict comple- 
mentarity into the framework of linear fractional optimization. We prove the 
existence of a strict complementary solution and show that all such solutions 
induce unique optimal partitions for the sets of indices of nonnegative vari- 
ables. To geometrically interpret the strict complementarity, we equivalently 
represent primal and dual optimal sets by two nonnegative polyhedral sets 
that are described only by equality constraints. Then we prove that each pair 
of relative interior points of these representing sets is a strict complementary 
solution, and vice versa. By this result, we deal with the problem of identify- 
ing a strict complementary solution. Specifically, we turn this problem to the 
equivalent problem of identifying a maximal element of a nonnegative poly- 
hedral set. Then, by applying the technique of finding a maximal element of 
a nonnegative polyhedral set, we develop two linear optimization approaches 
with different strategies for finding a strict complementary solution in linear 
fractional optimization. 

Our first approach identifies a strict complementary solution by solving 
two LPs and requires knowing the optimal objective value of the given LFP. 
In contrast, our second approach involves solving a single (but larger) LP 
and does not need the per-knowledge of the optimal objective value. As our 
proposed approaches are linear optimization based, they allow for applying 
the ordinary simplex algorithm of linear optimization to identify a strict com- 
plementary solution in linear fractional optimization. Nonetheless, regarding 
the preference on using the proposed approaches, note that each of the LPs 
developed in our first approach has less number of constraints than the LP 
of our second approach. Therefore, the use of our first approach is particu- 
larly recommended in situations, where only primal or dual part of a strict 
complementary solution needs to be found. For example, while the SBM 
model of Tone [26] is used for the measurement of efficiency in the field of 
data envelopment analysis, the global reference set of an inefficient decision 
making unit can be identified by either the primal or dual part of a strict 
complementary solution. However, solving the primal SBM model is recom- 
mended because the number of its constraints are mostly less than that of 
its dual. 

The approaches developed in our paper open up a number of further 
research avenues. First, it should be interesting to extend the results pro- 
posed in linear optimization literature on the use of strict complementarity 
for post-optimality analysis [5, 13] to linear fractional optimization. Second, 
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an interesting method for finding a strict complementary solution in linear 
optimization is to apply the so-called Balinski-Tucker tableau [2]. From our 
contribution, it is found that this method can be used to generate (indirectly) 
a strict complementary solution in linear fractional optimization. Therefore, 
it is worth exploring modification of the Balinski-Tucker method so that a 
strict complementary solution is directly obtained. 


Appendix A 


The following computer program written in GAMS identifies a strict com- 
plementary solution for the primal LFP (8) and its dual (9). Making this 
program applicable for any LFP in the general form of (1) just requires 
modifying “Sets”, “Table A(i,7)”, “Parameters”, “Alpha,” and “Beta” in 
Lines 1-23. 


1 Sets 

2 i row number of matrix A /i1*i2/ 
3 j column number of matrix A /j1*j2/; 
4 

5 Table A(i,j) 

6 ji j2 
7 il 2 1 

8 i2 -2 Ly 
9 

10 Parameters 

11 b/it 6 
12 i2 2/ 
13 c/ji 6 
14 j2 3/ 
15 d/ji 5 
16 j2 2/; 
17 

18 Scalars 

19 Alpha 

20 Beta; 

21 

22 Alpha=6; 

23 Beta =5; 

24 


25 File ProgSC / Results.txt /; 
26 Put Progsc; 


28 3 2 2 2k 2 2 2 2 2 2 a 2 2 2 2 2 2 2 2 2 2 2 ke 2 2 2 2 2k 2 2 i ak 2 i 2 i 2k 2c a 2 a 2k 2 a 
29 ¥*Stage 1: Solving program (3) 


30 

31 Free Variables 
32 Theta; 

33 

34 Positive Variables 
35 xbar (j) 

36 t 3 

37 

38 Scalar 

39 ThetaStar; 
40 


41 Equations 
42 Obj 
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43 Cont 

44 Con2; 

45 

46 Obj.. Theta =E= Sum(j, c(j)*xbar(j)) + Alpha*t; 

47 Coni(i).. Sum(j, a(i,j)*xbar(j)) =L= b(i)*t; 
48 Con2.. Sum(j, d(j)*xbar(j)) + Beta*t =E= 1; 
49 

50 Model MainLFP / Obj, Coni, Con2 /; 

51 

52 Put /'Finding the Optimal Obj. Value (ThetaStar)'; 

530 Put /'----- nnn nnn nnn nn nnn nnn nnn nnn nnn nana "Le 


54 Option LP=CONOPT; 
55 Solve MainLFP using LP Maximizing Theta; 


56 Put 'Obj = ':>6; Put Theta.L:<10:3; 

57 ThetaStar=Theta.L; 

88 Put /'--- oar nnn sn "4; 
59 


60 *End of Stage 1 
61 3 2 2 2 2 2 2 2 2 2 2 2 2 2 ke 2 2 ke 2k 2 ke 2 ke 2k 2 ik 2c a 2 ak 2 a 


62 
63 22 2 2 2 2 2 2 2 i 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 i 2 2 2 2 2 2 ie 2 2 2 ke 2 2 ie 2 2 2 ke 2 2K 2 2 2K 2 2k 2 OK 


64 *First approach: Solving program (18) 


65 

66 Positive Variables 
67 xbar1(j) 
68 xbar2(j) 
69 ubar1 (i) 
70 ubar2(i) 
71 y1(i) 

72 y2(i) 

73 v1(j) 

74 v2(j) 

75 wi 

76 w2 

77 P3 

78 

79 Free variable 
80 q3 

81 


82 xbar2.up(j) = 1; 
83 ubar2.up(i) = 1; 
84 y2.up(i) = 1; 

1 


85 v2.up(j) = 1; 
86 w2.up = 1; 

87 

88 Parameters 

89 XbarStar 
90 tStar 

91 UbarStar 
92 YStar (i) 
93 zStar 

94 vStar(j); 
95 

96 Equations 

97 ObjP 

98 ConP1 
99 ConP2 
100 ConP3 
101 ObjD 

102 ConD1 
103 ConD2 
104 ConD3; 
105 

106 ObjP.. Theta =E= Sum(j, xbar2(j)) + Sum(i, ubar2(i)) + w2; 
107 ConP1(i).. Sum(j, aCi,j)*(xbar1(j)+xbar2(j))) - b(i)*p 


+ ubari(i)+ubar2(i) =E= 0; 
108 ConP2.. Sum(j, d(j) ¥*(xbar1(j)+xbar2(j))) + Beta*p 
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- witw2 =E= 0; 

109 ConP3.. Sum(j, c(j) ¥*(xbar1(j)+xbar2(j))) + Alpha*p 
- (witw2)*ThetaStar =E= 0; 

110 

111 ObjD.. Theta =E= Sum(i, y2(i)) + Sum(j, v2(j)) + w2; 

112 ConD1(j).. Sum(i, a(i,j)*(y1(i)+y2(i))) + d(j)*q - viCj) 
- v2(j) - c(j)*(witw2) =E= 0; 

113 ConD2.. -Sum(i, b(i) ¥*(y1(i)+y2(i))) + Betatq 
- Alpha*(w1+w2) =E= 0; 

114 ConD3.. q - ThetaStar*(wi+w2) =E= 0; 

115 

116 Models Primal_SCSC / ObjP , ConP1, ConP2, ConP3 / 

117 Dual_SCSC / ObjD , ConD1, ConD2, ConD3 / ; 

118 


119 Solve Primal_SCSC using LP Maximizing Theta; 


120 XbarStar(j) = (xbar1.L(j)+xbar2.L(j))/(w1.L+w2.L); 
121 tStar = p.L/(w1.L+w2.L) ; 

122 UbarStar(i) = (ubar1.L(i)+ubar2.L(i))/(w1.L+w2.L); 
123 


124 Solve Dual_SCSC using LP Maximizing Theta; 


125 YStar(i) = (y1.L(i)t+y2.L(i))/(w1.L+w2.L); 

126 zStar = q.L/(w1l.L+w2.L) ; 

127 vStar(j) = (v1.L(j)+v2.L(j))/(w1.L+w2.L) ; 

128 

129 Put / /'Finding a SC Solution via Approach I'; 

130 Put /'----------------------------------------------- '/; 

131 Put ' Primal Dual '/; 

132 “Puts seer areas tS ers yt ae ee Oe A 

133  Loop(j, 

134 Put 'x(':>5; Put ord(j):<>3:0; Put ')= ':3; Put (XbarStar(j)/tStar) :<10:3; 
135 Put 'v(':>5; Put ord(j):<>3:0; Put ')= ':3; Put Vstar(j):<10:3; 
136 Put /; 

137 Js 

138 Put /; 

139 Loop(i, 

140 Put 'u(':>5; Put ord(i):<>3:0; Put ')= ':3; Put (UbarStar(i)/tStar) :<10:3; 
141 Put 'y(':>5; Put ord(i):<>3:0; Put ')= ':3; Put Ystar(i):<10:3; 
142 Put /; 

143 Ds 

144 Put '----------------------------------------------- Y11; 

145 


146 *End of First approach 
147 2 Fe 2 2 fe 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 i 2 2 2 2 ke 2 2 2 2 ke 2 2 2 2 2k 2 2k 2 2c a 2 ko 


148 
149 2 2 2 ke 2 2 2 2k 2 2k 2k 2 2 2 2k 2 2g 2k 2 2 2 2k 2 2k 2k 2 2 2 2k 2 2k 2k 2k 2 2 2k 2 2k 2k 2 2 2 2k 2 oi ke 2 2 2 2k 2 ok 2k 2 2 ok ke 2 ok 


150 *Second approach: Solving program (22) 


151 

152 Equations 

153 Obj PD 

154 ConPD ; 

155 

156 ObjPD.. Theta =E= Sum(j, xbar2(j)) + Sum(i, ubar2(i)) + Sum(i, y2(i)) 
+ Sum(j, v2(j)) + w2; 

157 ConPD.. Sum(j, c(j)*(xbar1(j)+xbar2(j))) + Alpha*p - q =E= 0; 

158 

159 Model PD_SCSC / ObjPD, ConP1, ConP2, ConD1, ConD2, ConPD/ ; 

160 

161 Solve PD_SCSC using LP Maximizing Theta; 

162 XbarStar(j) = (xbar1.L(j)+xbar2.L(j))/(w1.L+w2.L) ; 

163 UbarStar(i) = (ubari.L(i)+ubar2.L(i))/(w1.L+w2.L); 

164 Ystar(i) = (y1.L(i)t+y2.L(i))/(w1.L+w2.L); 

165 Vstar(j) = (v1.L(j)+v2.L(j))/(w1.L+w2.L) ; 

166 

167 Put / /'Finding a SC Solution via Approach II'; 

168. “Pith. f Saar cen stet aR eee Saat 4s 

169 Put ' Primal Dual '/; 


170 Put '----------------------------------------------- '/; 
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171 Loop(j, 

172 Put 'x(':>5; Put ord(j):<>3:0; Put ')= ':3; Put (XbarStar(j)/tStar) :<10:3; 

173 Put 'v(':>5; Put ord(j):<>3:0; Put ')= ':3; Put Vstar(j):<10:3; 

174 Put /; 

175 3 

176 Put /; 

177 Loop(i 

178 Put 'u(':>5; Put ord(i):<>3:0; Put ')= ':3; Put (UbarStar(i)/tStar) :<10:3; 

179 Put 'y(':>5; Put ord(i):<>3:0; Put ')= ':3; Put Ystar(i):<10:3; 

180 Put /; 

181 )3 

182 Put '----------------------------------------------- Yi] 

183 

184 »*End of Second approach 

185 2 2 2 2 2 i 22 2 2k 2 2 2k 2 2 2k 2 2 2 2k 2 2k 2k 2 2 2 2k 2 2k 2k 2 2 2 2k 2 2 2k 2k 2 i 2c 2 ok ke 2 2 ke 2k 2 ok ke 2k 2 ok ke 2 ok 
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